Benzene Adsorbed on Metals: Concerted Effect of Covalency and van der Waals 
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The adsorption of aromatic molecules on metal surfaces plays a key role in condensed matter 
physics and functional materials. Depending on the strength of the interaction between the molecule 
and the surface, the binding is typically classified as either physisorption or chemisorption. Van der 
Waals (vdW) interactions contribute significantly to the binding in physisorbed systems, but the 
role of the vdW energy in chemisorbed systems remains unclear. Here we study the interaction 
of benzene with the (111) surface of transition metals, ranging from weak adsorption (Ag and 
Au) to strong adsorption (Pt, Pd, Ir, and Rh). When vdW interactions are accurately accounted 
for, the barrier to adsorption predicted by standard density functional theory (DFT) calculations 
essentially vanishes, producing a metastable precursor state on Pt and Ir surfaces. Notably, vdW 
forces contribute more to the binding of covalently bonded benzene than they do when benzene is 
physisorbed. Comparison to experimental data demonstrates that some of the recently developed 
methods for including vdW interactions in DFT allow quantitative treatment of both weakly and 
strongly adsorbed aromatic molecules on metal surfaces, extending the already excellent performance 
found for gas-phase molecules. 
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I. INTRODUCTION 

The adsorption of aromatic molecules at transition- 
metal surfaces is important for fundamental and ap- 
plied surface science studies^— and these systems show 
promise as components in (opto)-electronic devices! 4 - In 
the case of weak overlap of electron orbitals between the 
adsorbate and the substrate surface, the ubiquitous van 
der Waals (vdW) interactions is frequently the only force 
that binds the molecule to the surface. This situation is 
typically referred to as physisorption. In the chemisorp- 
tion case, the covalent or ionic bonding dominates and 
the effect of vdW interactions on the overall strength of 
adsorption is typically assumed to be weak. In this study, 
we challenge this conventional view, by demonstrating 
the significantly larger contribution of vdW energy to 
the stabilization of strongly adsorbed benzene on (111) 
surfaces of Pt, Pd, Rh, and Ir metals when compared to 
physisorption on Ag(lll) and Au(lll) surfaces. 

Whereas the role of vdW forces in the binding of atoms 
and molecules in the gas phase is reasonably well under- 
stood, at solid surfaces our understanding remains far 
from complete. Indeed, until recent developments (see, 
e.g., Refsja— ) for efficiently incorporating the long-range 
vdW energy within density functional theory (DFT) cal- 
culations it was not possible to determine the role of the 
vdW energy for extended systems and adsorption pro- 
cesses A large majority of previous theoretical work 
on vdW interactions mainly focused on weakly bound 



systems. 



Typical examples include benzene (Bz) ad- 



sorbed on the Ag(lll) and Au(lll) surfaces,— — and no- 
ble gases on the Cu(lll), Ag(lll), Pt(lll), and Pd(lll) 
surfaces. 10 ' 19 " — A unifying aspect of these studies is the 



observation that the inclusion of vdW interactions into 
standard DFT within the generalized gradient approx- 
imation (GGA) often brings a large increase in bind- 
ing, and results in a much better agreement with ex- 
perimental adsorption distances and energies. However, 
the vdW forces can also have a qualitative impact on 
the adsorption process. One particularly interesting ex- 
ample was reported by Bliigel's group, showing that the 
vdW forces are the key ingredient to trigger the bind- 
ing of pyridine on Cu(llO) from physisorption to weak 
chemisorptioni22 Mittendorfer et al£2- reported a novel 
mechanism for graphene adsorption on Ni(lll), where 
weak covalent and vdW interactions lead to two differ- 
ent minima in the binding curve. Similar results were 
found by Li et al.2^ Another example was shown in our 
recent work on the isophorone molecule (C9H14O) at the 
Pd(lll) surface, which illustrated that the binding struc- 
ture and the dehydrogenation pathway in this system 
can be predicted only after accounting for vdW interac- 
tions^ The vdW interactions were also shown to play a 
role in the chemisorption of benzene on the Si(100) sur- 
facei 26 ' 27 However, in this case the vdW-DF method leads 
to a smaller adsorption energy than the pure PBE func- 
tional^ while the PBE+vdW method predicts a vdW 
contribution of 0.5 eV<2I 

In this work, we demonstrate the significant concerted 
effect of covalent bonding and vdW interactions for ben- 
zene interacting with metal surfaces, leading to qualita- 
tive changes in the adsorption behavior when vdW inter- 
actions are accurately treated (see Figure 2). In particu- 
lar, our calculations predict a metastable precursor state 
for benzene on Pt(lll) in agreement with the experimen- 
tal findings, 28 along with a peculiar "phase transition" 
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behavior of the projected HOMO/LUMO occupations 
of the benzene molecule. Comparison to experimental 
data demonstrates that recently developed methods for 
including vdW interactions in DFT— ii2 allow quantitative 
treatment of both weakly and strongly adsorbed aromatic 
molecules on metal surfaces, extending the already excel- 
lent performance found for gas-phase molecules 



II. COMPUTATIONAL METHOD 

We used two different vdW-inclusive approaches in 
the present work: a newly developed PBE+vdW s 
method, 10 as implemented in the FHI-aims all-electron 
code; 30 and the optB88-vdW method. 8 as implemented 
in the VASP code^i^ The PBE+vdW surf approach 
includes screened vdW interactions (beyond the pair- 
wise atom-atom approximation) to study adsorbates 
on surfaces, by a synergetic linkage of the PBE+vdW 
method 6 for intermolecular vdW interactions with the 
Lifshitz-Zaremba-Kohn theory 3 - 3 - for the dielectric screen- 
ing within the metal substrate. While the PBE+vdW surf 
approach leads to accurate results in the asymptotic limit 
by construction, it uses a short-range damping function 
with one adjusted parameter. The optB88-vdW method 
is a modified version of the vdW-DF functional^ by us- 
ing an empirically optimized optB88-like exchange func- 
tional. Both PBE+vdW surf and optB88-vdW methods 
can accurately describe intermolecular interactions with 
mean absolute relative errors on the order of 9%&22 com- 
pared to coupled-cluster dimer binding energies for the 
S22 molecular database. Less is known about the per- 
formance of these methods for solids and weakly ad- 
sorbed molecules on surfaces, although encouraging re- 
sults have been reported for a few condensed matter sys- 
tems jioiLHiM j g i m p 0r t a nt to benchmark the newly 
developed methods on a wider set of condensed mat- 
ter systems, especially because the PBE+vdW surf and 
optB88-vdW methods are based on very different approx- 
imations. For comparison purposes, calculations using 
the vdW-DF functional and its second version (vdW- 
DF2)i£ were gj so carr i e(1 ou t for the Bz/Pt(lll) and 
Bz/Au(lll) systems. 

The FHI-aims code was employed for the 
PBE+vdW surf , PBE+vdW, PBE, 36 and local-density 
approximation (LDA) calculations. We used the "tight" 
settings, including the standard numerical atom-centered 
orbitals (NAO) basis set "tier2" for H and C, and "tierl" 
for transition metals. For all structural relaxations, 
we used a convergence criterion of 10~ 2 eVA^ 1 for the 
maximum final force. Also a convergence criteria of 10 -5 
electrons for the electron density and 10 -4 eV for the 
total energy of the system were utilized for all compu- 
tations. The scaled zeroth-order regular approximation 
(ZORA) was applied for treating relativistic effects. 
Using these settings, the accuracy in determining the 
binding energy and equilibrium distance is better than 
0.01 eV and 0.005 A, respectively. 



For bulk lattice constant calculations, we used a 
Monkhorst-Pack 3 -! grid of 16 x 16 x 16 fc-points. The 
lattice constants of the bulk metals have been obtained 
by using the Birch-Murnaghan equation of state fit to 
DFT cohesive energy curves^ Using the respective lat- 
tice constants from each method, we built up a 6-layer 
slabs with a (3 x 3) unit cell, with no reconstruction of 
Pt(lll) and Au(lll). Each slab was separated by a 20 
A vacuum. The vdW interactions between metal atoms 
were also considered when performing the relaxations. 
We constrained the bottom four metal layers while fully 
relaxed the molecule and the uppermost two metal lay- 
ers during geometry relaxations. For slab calculations, 
we used a 6 x 6 x 1 fc-points mesh. 

For the calculation of the binding curves, we changed 
the adsorption height d of Bz, which is evaluated relative 
to the position of the unrelaxed topmost metal layer. For 
each structure, we fixed the z coordinates of the carbon 
backbone and the metal atoms in the bottommost four 
of the employed six layer surface model. 

The VASP code was employed for the optB88-vdW, 
vdW-DF, and vdW-DF2 calculations. Inner electrons re- 
placed by the projector augmented wave (PAW) method, 
whilst the monoelectronic valence electrons were ex- 
panded in plane- waves with a E cnt - R = 500 eV. For slab 
calculations, we used a 4 x 4 x 1 fc-points mesh. For metal 
supercell, we used a (3 x 3) unit cell with 6 atomic layers 
(3 bottom layers fixed to the corresponding bulk optimal 
position for each method). Dipole correction was applied 
along the direction perpendicular to the metal surface. 
Geometry optimizations were performed with a residual 
force threshold of 0.03 eVA" 1 . 



III. RESULTS 

The typical strongly bound Bz/Pt(lll) system (ad- 
sorption energy 1.57-1.91 eV— ) and the typical weakly 
bound Bz/Au(lll) system (adsorption energy 0.73-0.87 
eV^) are used first to demonstrate our point. Accurate 
experimental data is available for both of these systems, 
enabling direct quantitative verification of our theoreti- 
cal calculations. To demonstrate the differences in the 
adsorption mechanism, we explore the potential-energy 
surface (PES) for Bz on the Pt(lll) and Au(lll) sur- 
faces. We place a single Bz molecule at the eight high- 
symmetry adsorption sites of the (111) metal surface^ - 
followed by geometry relaxation. The adsorption geome- 
tries and energies for Bz on Au(lll) and Pt(lll) at the 
preferable adsorption site are shown in Fig. 1 and Table 
I. Already here one can clearly distinguish the different 
nature of bonding for the adsorption of Bz on Pt(lll) 
and Au(lll). Irrespective of the functional used (PBE, 
PBE+vdW, and PBE+vdW surf ), the bri30° is the most 
preferable site for Bz/Pt(lll), with an angle of 30° be- 
tween the C-C and Pt-Pt bonds, see Fig. 1. This re- 
sult is consistent with previous periodic slab GGA cal- 
culations as well as low-energy electron diffraction 
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(LEED)^L and scanning tunneling microscopy (STM) 46 
experiments. Moreover, the PES shows a corrugation of 
1.33 eV for Bz/Pt(lll) when using PBE+vdW surf . In 
contrast, the PES for Bz/Au(lll) is found to be flat, 
with only 0.04 eV corrugation. This result further jus- 
tifies the STM observations that even at a temperature 
of 4 K, Bz molecules are capable of diffusing over the 
Au(lll) terraces4i 
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FIG. 1. Adsorption structures of the Bz/Pt(lll) system and 
Bz/Au(lll) system, both at the so-called bri30° adsorption 
site (see text). We carried out extended periodic calculations, 
but only a small part of the supercell is shown. Six metal lay- 
ers were used but only the topmost three layers are depicted in 
the figure. The indicated distances (A) are obtained based on 
the PBE+vdW surf optimized structures. Gray, yellow, cyan, 
and white spheres represent Pt, Au, C, and H atoms, re- 
spectively. Optimized lattice constants were used for every 
method. 

The analysis of the equilibrium distances and ad- 
sorption energies in Table I demonstrates that both 
PBE+vdW surf and optB88-vdW methods lead to an ex- 
cellent agreement with the available experimental data. 
For the Bz/Pt(lll) system, the PBE+vdW surf adsorp- 
tion energy of 1.96 eV is close to that from optB88-vdW 
(1.84 eV) and both methods agree with the measured 
calorimetry values at 0.7 ML (1.57-1.91 eV, the same cov- 
erage used for DFT calculations)^ The PBE+vdW surf 
adsorption energy converges to 2.18 eV with increasing 
surface cell size, within the error bar of calorimetry mea- 
surements in the limit of zero coverage (1.84-2.25 eV)<2£ 
Note that the exclusion of the vdW interactions in the 
strongly adsorbed Bz/Pt(lll) system would lead to a sig- 
nificant reduction in the binding energy (0.81 eV from 
PBE), in disagreement with the experimental data. The 
adsorption energies computed using the vdW-DF and 
vdW-DF2 methods are even smaller than those calcu- 
lated with PBE. The adsorption energy for Bz adsorbed 



TABLE I. Comparison of adsorption energy (-E a d) and av- 
erage perpendicular heights (rfcM and (Zhm for carbon-metal 
and hydrogen- metal, respectively) between DFT calculations 
and experimental data for Bz on Pt(lll) and Au(lll). The 
distances are referenced to the average positions of the relaxed 
topmost metal atoms. The adsorption energy E^d is defined 
as i? a d = -(EAdSys-EMe-E-Bz), where the subscripts AdSys, 
Me, and Bz denote the adsorption system, the clean metal 
substrate, and the isolated Bz molecule, respectively. 



System 


Method 


-Bad feVl 


dcM Al 


dfJ'M [Al 




PBE+vdW surt 


1.96 


2.08 


2.51 




optB88-vdW 


1.84 


2.12 


2.53 


Bz/Pt(lll) 


vdW-DF 
vdW-DF2 


0.77 
0.34 


2.16 
2.20 


2.57 
2.65 




PBE 


0.81 


2.10 


2.54 




LDA 


2.30 


2.05 


2.47 




Experiment 


1.57-1.91 a 


2.02±0.02^ 






PBE+vdW surt 


0.74 


3.05 


3.04 




optB88-vdW 


0.79 


3.23 


3.23 


Bz/Au(lll) 


vdW-DF 


0.59 


3.44 


3.42 


vdW-DF2 


0.56 


3.29 


3.27 




PBE 


0.15 


3.62 


3.62 




LDA 


0.49 


2.83 


2.82 




Experiment 


0.73-0.87 c 


2.95-3.10 d 





a Heat of adsorption measured with calorimetry, at the same 
coverage (0.7 ML) used for the DFT calculations^ The error 
estimates of ±10% are taken from the referenced Recent work 
suggests reduced errors of ±5%— 

b LEED experiment — 

c TPD experimental^ 

d Deduced data based on the experimental workfunction for Bz 
on Au(lll) and adsorption distance for pentacene on 
Au(lll)i^^ 



on the Au(lll) surface is considerably smaller than that 
of Bz/Pt(lll). Also for Bz/Au(lll), the PBE+vdW surf 
adsorption energy (0.74 eV) agrees very well with both 
the optB88-vdW result (0.79 eV) and the experimental 
temperature-programmed desorption (TPD) data at 0.1 
ML (0.73-0.87 eV)^®. We conclude that PBE+vdW surf 
and optB88-vdW methods yield quantitative agreement 
with experimental adsorption distances and energies for 
both weakly and strongly adsorbed Bz molecule. In con- 
trast, LDA calculations are not systematic, underbinding 
for Au(lll) and overbinding for Pt(lll). 

Deeper insight into the mechanism of Bz adsorption 
can be gained upon analysis of the binding energy curves, 
End(d), in Fig. 2. The binding energy curves exhibit 
several characteristic effects. With decreasing distance 
the binding energy of the adsorbate system increases, 
determined mainly by vdW interactions, and here (for 
rf>3.5 A) Au(lll) and Pt(lll) show very similar behav- 
ior. In both cases the calculations show a small broad- 
ening of the energy levels. The fully occupied rf-band 
of Au is obviously stiffer than the partially empty ri- 
band of Pt. In fact, for the latter the Pauli repulsion can 
be weakened by the rearrangement of rf-electron density 
(a similar effect has been investigated in previous liter- 
ature^). As a consequence, the Bz molecule gets closer 
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FIG. 2. Top: Adsorption energy -E a d as a function of the adsorption height d for Bz on Pt(lll) (a) and on Au(lll) (b) 
from the PBE and PBE+vdW smf methods (the carbon backbone height d from the surface is kept fixed). The experimental 
binding distances and adsorption energies are indicated by yellow intervals. Bottom: Integrated projected molecular density of 
states^ for the HOMO-1, HOMO, LUMO, and LUMO+1 orbitals of the benzene molecule as a function of d for Bz on Pt(lll) 
(c) and on Au(lll) (d). The inset in panel (c) shows a side view of the electron density difference, which was obtained by 
subtracting electron density of isolated molecule and clean surface from an electron density plot of the entire adsorbed system, 
upon Bz adsorption on Pt(lll) at d=2.08 A (red = electron depletion, blue = electron accumulation). For the same value of 
the isosurface (0.04 A -3 ), the electron density difference for Bz/Au(lll) at d=2.08 A is significantly weaker, see the inset in 
panel (d). 



to the surface of Pt and the HOMO and LUMO levels 
of the combined system broaden and hybridize notice- 
ably. This goes together with significant electron trans- 
fer: The HOMO and HOMO-1 orbitals of Bz molecule 
get partially depleted and the LUMO and LUMO+1 or- 
bitals become partially filled. This behavior (broadening, 
shift, hybridization of levels, and electron transfer) is a 
clear signature of the covalent interaction for Bz/Pt(lll). 
Thus, at the adsorption geometry the wave-function has 
attained a qualitatively new character. Figure 2 shows 
that this character change sets in for Pt at a distance of 
3.1 A. At 2.6 A nearly a full electron has been transferred 
from the HOMO and HOMO-1 levels to the LUMO and 
LUMO+1 levels, and in the total energy we observe a 
"phase transition behavior" (cf. the peak at 2.6 A). Fi- 
nally, at the equilibrium geometry the electron transfer 
(rearrangement) is as large as ~1.1 electrons. For Au sur- 
face the process is much weaker and - not surprisingly 
- a covalent contribution to the adsorption process re- 
mains negligible. Thus, the vdW attraction governs the 
interaction. 

Further inspection of the electron density difference 
at the strongly bound minimum for Bz/Pt(lll) in Fig. 



2(c) demonstrates the rather strong hybridization be- 
tween the HOMO/LUMO orbitals of Bz and the d z 2 or- 
bitals of the Pt(lll) atoms. For the same adsorption 
height, the electron density difference for Bz/Au(lll) is 
weak (see Fig. 2(d), inset). The presence of two minima 
for Bz/Pt(lll) resembles the recently studied bonding 
of graphene on Ni(lll)i 23 ' 24 However, the adsorption of 
Bz on Pt(lll) exhibits a different feature. In fact, Bz 
is exothermically bound on Pt(lll) already when using 
PBE without vdW interactions, while the PBE adsorp- 
tion energy is endothermic for graphene on Ni(lll). Evi- 
dently, the functionalization of aromatic molecules would 
allow to control the position and stability of the two ad- 
sorption minima on metallic surfaces. 

Interestingly, while covalency is crucial for the 
Bz/Pt(lll) bonding character, energetically the vdW 
contribution is in fact significant. Upon inclusion of 
vdW interactions, the binding behavior is strongly mod- 
ified - the barrier to adsorption vanishes, and a pre- 
cursor physisorption state emerges for Bz/Pt(lll). The 
PBE+vdW surf method lowers the adsorption energy from 
0.50 eV (pure PBE value) to 1.65 eV in Fig. 2. Thus the 
final adsorption results from a strongly concerted, syn- 
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TABLE II. Adsorption energies i? a d (eV) of Bz adsorbed on 
(111) surfaces of Ag, Pd, Rh, and Ir. 



System 


PBE 


PBE+vdW surt 


optB88-vdW 


Bz/Ag(lll) 


0.09 


0.75 


0.72 


Bz/Pd(lll) 


1.17 


2.14 


1.91 


Bz/Rh(lll) 


1.48 


2.52 


2.27 


Bz/Ir(lll) 


1.10 


2.24 


2.09 



ergistic effort. Upon comparing the binding curves for 
Bz/Pt(lll) and Bz/Au(lll) we see that the vdW con- 
tribution (due to vdW surf ) for Bz/Pt(lll), 1.15 eV, is 
even stronger that that for Bz/Au(lll), 0.68 eV. The 
screened Bz/surface C3 vdW coefficient is essentially the 
same for Pt(lll) and Au(lll) surfaces (2.17 and 2.02 
hartree bohr 3 , respectively). Therefore, we conclude that 
the larger contribution of the vdW energy in the case of 
covalent bonding comes from the rather short adsorption 
distance of the Bz molecule from the surface. 

Our conclusions hold in general for the adsorption of 
Bz on other transition metal surfaces. For Bz/Ir(lll), 
the binding curve shows the same characteristic features 
as for Bz/Pt(lll) in Fig. 2. For Bz adsorbed on the 
Pd(lll), Rh(lll), and Ir(lll) surfaces, the vdW en- 
ergy contributions from the PBE+vdW s method are in 
the range of 0.97-1.21 eV, greater than those for Bz ph- 
ysisorbed on Ag(lll) and Au(lll) (0.68-0.82 eV). Even 
larger vdW energies are found in more complex pol- 
yaromatic adsorption systems. For instance, the vdW 
energy is determined to be 1.77 eV for naphthalene 
(CioH 8 ) on the Pt(lll) surface with (5 x 4) unit cell. 
Also for this case the calculated adsorption energy from 
PBE+vdW surf (2.91 eV) is within the experimental er- 
ror bars (2.80-3.42 eV)^ For anthracene (C u Hi ) on 



the Pt(lll) surface with (6 x 4) unit cell, the adsorption 
energy contributed by vdW interactions (2.42 eV) largely 
exceeds that determined from the PBE functional (1.38 
eV). 



IV. CONCLUSION 

In summary, we have demonstrated that the inclusion 
of vdW interactions qualitatively changes the adsorp- 
tion behavior for benzene strongly interacting with (111) 
metal surfaces. The vdW energy in Bz/Pt(lll), a typical 
strongly adsorbed system, is almost 0.5 eV greater than 
that in Bz/Au(lll), a typical physisorbed system. The 
bonding mechanism of Bz/Pt(lll) stems from a syner- 
gistic effort of covalent bonding and vdW interactions, 
and it is characterized by a peculiar "phase transition" 
behavior in the projected HOMO/LUMO occupations of 
the benzene molecule. Our findings for Bz adsorbed on 
Ag, Au, Pt, Pd, Rh, and Ir surfaces indicate that DFT 
calculations with dispersion interactions are essential for 
both weakly and strongly bound molecules on surfaces. 
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